Supercurrent-carrying density of states in diffusive mesoscopic Josephson weak links 
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Recent experiments have demonstrated the nonequilibrium control of the supercurrent through 
diffusive phase-coherent normal-metal weak links. The experimental results have been accurately 
described by the quasiclassical Green's function technique in the Keldysh formalism. Taking into 
account the geometry of the structure, different energy scales and the nonidealities at the interfaces 
allows us to obtain a quantitative agreement between the theory and the experimental results in 
both the amplitude and the phase dependence of the supercurrent, with no or very few fitting 
parameters. Here we discuss the most important factors involved with such comparisons: the ratio 
between the superconducting order parameter and the Thouless energy of the junction, the effect 
of additional wires on the weak link, and the effects due to imperfections, most notably due to the 
nonideal interfaces. 

PACS numbers: 74.50.+r, 73.23.-b, 74.40.+k, 74.80.Fp 



I. INTRODUCTION 

Many quantum phenomena in many-body systems are 
based on probing the spectrum of states corresponding to 
the desired observable, the states being filled according to 
an appropriate distribution function. A similar viewpoint 
can be taken also on the Josephson effect: supercurrent 
is carried by states in the weak link and their occupation 
is determined by a distribution function antisymmetric 
between the electron and hole spaces. This aspect is 
directly reflected in the mathematical structure of the 
supercurrent formula derived from the Keldysh Grecn's- 
functions method. 1-3 Such an approach has been taken 
in some recent experiments 4-9 controlling the Josephson 
effect in phase-coherent normal-metal wires through the 
control of the distribution function by an injection of nor- 
mal quasiparticle current. One of the most remarkable 
results of these experiments is the inversion of the sign of 
the supercurrent for a given phase difference across the 
weak link when the junctions turn into a 7r state. 

Quantitative fit to the experimentally obtained re- 
sults has been very successful for the equilibrium 
supercurrent 10 using the equilibrium quasiclassical the- 
ory. In the nonequilibrium case, detailed knowledge 
of the relaxation mechanisms controlling the shape 
of the interactions, but also the precise spectrum of 
supercurrent-carrying states, is required. 6 ' 8 Previously 1 
for the calculation of this spectrum, one has assumed a 
two-probe setup with some idealized conditions on the 
length scales and on the nature of the interfaces. In this 
paper, we systematically investigate the spectrum of this 
current-carrying density of states, or spectral supercur- 
rent, show how it is calculated, and how it depends on 
the length of the weak link, presence of additional termi- 
nals, or on the nonidealities in the interfaces between the 
normal-metal weak link and the superconductors. We 
also discuss the current-phase relation of such a system: 
at low temperatures, it can be far from sinusoidal, and at 



certain conditions, its period can even be halved. 11 We 
focus on the diffusive limit where the dimensions d of the 
weak link arc much greater than the elastic mean free 
path I. This is the typical limit for most normal- metal 
weak links. The corresponding ballistic limit d <C I has 
been extensively described in the literature 12-18 in terms 
of Andreev bound states (ABS). We show qualitatively a 
connection between the discrete ABS and the continuous 
diffusive-limit spectral supercurrent. 

This paper is organized as follows. After this intro- 
duction, Sec. II introduces to the theoretical formalism 
which is based on the real-time Usadel equation for the 
quasiclassical Green's function in the diffusive limit. 19 ' 20 
In the case of nonideal interfaces or in multiterminal ge- 
ometries, the boundary conditions to these functions are 
also essential. Understanding the results of the follow- 
ing sections does not require a detailed reading of this 
part but it is enough to grasp the idea of the relation 
of the spectral supercurrent and the observable one. In 
Sec. Ill we look how the spectral supercurrent depends 
on the length of the weak link compared to the super- 
conducting coherence length and separate two extreme 
cases. In the limit of a short junction where the co- 
herence length is much longer than the weak link, one 
obtains an analytical solution for the spectral supercur- 
rent without further approximations. The current-phase 
relation in diffusive normal-metal weak links is consid- 
ered in Sec. IV. We show how, especially at low temper- 
atures, higher harmonics appear in addition to the usual 
sinusoidal phase dependence and indicate how the period 
can be halved in a nonequilibrium situation. Section V 
considers the effect of additional normal-metal terminals 
on the current-carrying density of states, and in Sec. VI, 
we discuss how nonidealities in the normal metal — su- 
perconductor (NS) interface change its shape. Finally, in 
Sec. VII, we summarize the main results. 

To be specific, we consider the structure shown in Fig. 
1. The main wire with length L and cross section A w be- 
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we adopt the Keldysh real-time formalism 21 and hence 



FIG. 1: Multiterminal SNS Josephson junction. The weak 
link consists of a phase coherent normal-metal wire of length 
L, cross section A w , and normal-state conductivity a w , yield- 
ing a normal-state resistance Rm — L/a w A w . Additional 
normal- metal wires of length L c , total cross section A c , 
and normal-state conductivity a c , i.e., with resistance R c — 
L c /a w A w , called the control wires, are connected to the cen- 
ter of the weak link, and from their other end, to normal 
reservoirs. 



tween the superconductors forms the weak link whereas 
the additional wires with length L c and integrated cross 
section A c are used for the control of the distribution 
functions and therefore referred to as the control wires. 
We assume that the superconducting and normal reser- 
voirs are much larger than the weak link and the control 
wire, such that the Green's functions describing them 
take their bulk values very close to the interfaces. Fur- 
thermore, we assume that the width of the control wires 
is much smaller than the length L of the weak links. This 
allows us to consider the wires as quasi-one-dimcnsional 
structures by assuming translational invariance in the 
transverse directions. 



II. THEORETICAL BACKGROUND 

Circuits composed of normal and superconducting 
metals in the diffusive limit (dimensions larger than the 
elastic mean free path I) are effectively described in terms 
of the quasiclassical Green's functions G satisfying the 
Usadcl equations 19 ' 20 



DV(GVG)= -i(E + iT)t 3 + A,G 



(1) 



where D = ^vpl is the diffusion constant, E is energy 
relative to the chemical potential of the superconductors 
(which is assumed to be the same for all S terminals), 
T describes a small inelastic scattering rate, and A the 
superconducting pair potential (we set Ti = 1 through- 
out). Since we aim to describe nonequilibrium effects, 
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All of the submatrices denoted by a hat (G R , etc.) are 
2x2 matrices in Nambu particle-hole space, in particular, 
T3 is the third Pauli matrix and A has the form 



A(x) 
A*(x) 



(3) 



The pair potential A(x) can in principle be obtained from 
a self-consistency relation. 20,22 However, since we con- 
sider only superconducting reservoirs much wider than 
the weak link, we adopt the usual step-function form for 
A(x) (finite constant in the superconductors, zero in the 
normal- metal wires). 

In addition to Eq. (1), Usadel Green's function sat- 
isfies a normalization condition (5 =1. Therefore 
it can be parametrized with four scalar parameters as 
follows. 20 The Keldysh Green's function G K describ- 
ing the occupation numbers of different quantum states, 
i.e., the (non)equilibrium state of the system can be ex- 
pressed with two real distribution functions /x and 
as G K = G R {f L + f T T 3 ) - (f L + f T T 3 )G A whereas the 
retarded and advanced Green's functions, G R and G A , 
describing the spectral properties which do not directly 
depend on the distribution functions are 



G 



cosh(6*) sinh(6>) exp(i%) 

sinh(#) cxp(— i\) — cosh(6*) 



(4) 



and G A = -t 3 (G r )H 3 . Here 6{x;E) and x(x;E) are in 
general complex scalar functions. 

In what follows, we describe a quasi-one-dimensional 
situation, where the functions are assumed to vary only 
in one dimension x. Expressing the coordinate x in terms 
of the separation L of the superconductors between which 
the supercurrent flows, x = x'L, the spectral equations 
for G R ( A ^ read in a normal metal (A = 0) 

d 2 x ,6 = -2i(E' + iT') sinh(0) + i(cVx) 2 sinh(20), (5) 



j E = - smh 2 {e)d x >x, d x ,j E = 0. 



(6) 



Here, the prime over the (dimcnsionlcss) quantities de- 
notes the fact that the energies are expressed in the units 
of the Thouless energy Et = D/L 2 corresponding to the 
length L. Below, we tacitly assume all lengths and ener- 
gies expressed in these natural units even if not marked 
by a prime. The kinetic equations satisfied by the dis- 
tribution functions /l and fr are described, e.g., in Ref. 
20, where the part of the distribution function which is 
symmetric about the chemical potential of the supercon- 
ductors corresponds to Jt and the antisymmetric part to 
fi. These two components acquire different space and 
energy dependent diffusion coefficients due to the super- 
conducting proximity effect. 



If the interfaces to the reservoirs are ideal metallic, 
the parameters are continuous at the boundaries to the 
reservoirs and can be identified with the bulk values, 
8s = artanh(A/i?) and 9 m = in the superconduct- 
ing and normal-metal reservoirs, respectively. In general, 
e.g., if a supercurrent is driven through the system, there 
can be a phase difference, which we choose to be ap- 
plied symmetrically between the superconductors, such 
that in the left superconductor x — <fr/2 and in the right 
X = —<t>/2- Below, if not mentioned otherwise, we choose 
4> = n/2, which typically yields a supercurrent close to 
the critical current of the junction. 

Nonideal interfaces with reduced transmissivities are 
not directly described by the Usadel equation, because 
they are of microscopic, atomic-scale thickness. They 
can, however, be taken into account using boundary con- 
ditions derived by Zaitsev 23 for Eilenberger Green's func- 
tions (valid independent of the mean free path) and later 
simplified in the diffusive limit by Kuprianov and Lu- 
kichev for a tunneling case 24 and Nazarov for a general 
interface, 25 described by a scattering matrix. For an in- 
terface characterized by the transmission eigenvalues T n , 
the Green's functions G\ on the right-hand side and C? 2 
on the left-hand side of the interface satisfy 25 ' 26 



]_ w 2 w w 

a N AiGid x Gi = a N A 2 G 2 d x G 2 



N 

2e 2 



E 



T n [Gi, G2] 



(7) 



tt ^ 4 + T„({Gi,G 2 }-2)' 



evaluated at the position of the interface. In most cases, 
the individual transmission eigenvalues are not known, 
but since typical interfaces contain a huge number of 
channels, it is enough to integrate over the probabil- 
ity distribution of the eigenvalues to obtain the desired 
boundary condition. 

In the case of a tunneling interface (where all the 
transmission eigenvalues of the interface are small) the 
boundary conditions between the parametrized functions 
in wires 1 and 2 reduce to 20,24 

d x 9i = [sinh(6*i) cosh(0 2 ) — sinh(6> 2 ) cosh(6>i) cos(Ax)]/n,, 



sinh 2 (6>i)<9 x xi = sinh(0i) sinh(# 2 ) sin(Ax)/r 6 . 



(8) 
(9) 



Here, Ax = Xi ~ X2 and 6> 1(2) = 9(x^ ( ') and xi(2) = 

x( x t^) are tnc parameters 9 and x at the interface, 
x = Xb, but on the side of the wire 1 (2). The non- 
ideality of the interface is characterized by the ratio of 
its resistance Ri and of the weak-link resistance Rn, 
rj, = Ri / Rn and the derivatives point towards the wire 
1. In the case of a dirty interface, where the boundary 
condition is evaluated using the distribution function of 
the transmission eigenvalues corresponding to an inter- 
face with a random array of scatterers in a 2D layer, 27 



we get 



d x 6 x 



V2[sinh(6»i) cosh(6» 2 ) - sinh(6> 2 ) cosh(0i) cos(Ax)] 



sinh 2 (9i)d x xi 



r b V 



V2sinh(6> 2 ) sinh(0i) sin(Ax) 



(10) 
(11) 



Here we denoted the denominator V = 
y/1 + cosh(fl 2 ) cosh(fli) - sinh(0 2 ) sinh(#i) cos(Ax). 
This denominator reflects the contribution of open 
conduction channels which are not present in Eq. (9). 

Note that both types of boundary conditions indicate 
a form of a conservation of a spectral current over the 
interface, the second equation being the conservation of 
the spectral supercurrent je- 

In geometries with more than two terminals, we 
assume that narrow quasi-one-dimensional wires con- 
nect to each other at some point of the structure. 
Therefore we need to impose appropriate matching 
conditions. 20,25 ' 28,29 In this case, they are the continu- 
ity of the functions 9 and x and the conservation of the 
spectral currents. Assuming that the derivatives in the 
TV wires i = 1,...N with cross sections Ai and normal- 
state conductivities a l N point towards the crossing point 
at x c , we get 



e i { Xc ) = e j (x c ) vi,j = i,..jv, 

Xi{x c ) = Xj(x c ) = 1, -N, 

N 

Y^Ai(j^d x 6i{x c ) = 0, 



i=i 

N 



^AiG % N d x Xi{x c ) = 0. 



(12) 
(13) 

(14) 
(15) 



1=1 



In the last condition we used the continuity of the pa- 
rameters 9 across the crossing point. 

Below, we assume the system depicted in Fig. 1: two 
superconductors connected by "horizontal" mesoscopic 
normal wires to which we connect normal reservoirs by 
the "vertical" mesoscopic normal wires (labeling of the 
wires as in Fig. 1). When considering the supercurrent 
between the two superconductors, for the spectral equa- 
tions it is enough to treat any number of "vertical" wires 
by a single wire for which the product of a^A is simply 
the sum of these products in the individual wires. In the 
case that the dependence on the length L c of these wires 
becomes important, the smallest of them characterizes 
the situation the best. In this case, since there can be no 
supercurrent flowing to the normal reservoirs, Eq. (15) 
reduces to j E = Furthermore, for simplicity, we 

assume the system left-right symmetric, such that the 
part of the weak link in the left-hand side of the cross is 
similar to that in the right-hand side. 

Finally, the observable supercurrent is obtained from 
the solutions to the spectral and kinetic equations by 



En 



2eR 



N 



f 

J — ( 



dE'lm\j E {E')]f L {E'). (16) 
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In the reservoirs with voltage V with respect to the po- 
tential of the superconductors (which are assumed equal 
for both superconductors in order to avoid the ac Joseph- 
son effect), Jl obtains the form 



f L {E;V,T) = 



1 



tanh 



E + eV 
2k B T 



+ tanh 



E-eV 
2k B T 



(17) 

It can be shown 1,6 that, in the absence of inelastic inter- 
actions and for energies E < A, Jl remains constant 
throughout the control wires, and hence the reservoir 
value can directly be used for the calculation of the su- 
percurrent. 

In this paper, we will consider two limits for fa. 
These are the equilibrium finite-temperature limit, where 
/l = tanh(i?/2fcsT), and the zero-temperature nonequi- 
librium case when fp is driven in a normal-metal wire, 1 
f L = ${E - eV) - i?(-eV - E), where d(E) is the Heav- 
iside step function. 

The spectrum of supercurrent-carrying states typically 
consists of both the states carrying the supercurrent par- 
allel to the phase gradient and those carrying it in the 
opposite direction, 1,12 depending on their energy. Hence, 
by controlling the occupation of these states by the above 
stcplike distribution function, one is able to vary the sign 
of the observable supercurrent and, e.g., obtain the tt 
state. 

The form of the spectrum can be qualitatively under- 
stood by considering a ballistic (scattering-free) weak 
link. There, the quasiparticles form bound states 12-18 
which contain an Andreev reflection 30 at both NS inter- 
faces. Since the first reflection at the left is from hole- 
to particle (particle- to hole) -like states and the second 
at the right interface from particle- to hole (from a hole- 
to particle) -like states, the net result is a transfer of a 
Cooper pair from the left superconductor to the right 
(from right to left). Bound-state energies are found by 
requiring that the total phase the quasiparticles acquire 
within a single cycle is a multiple of 2n. This leads to 
(for E m < A) 



Et = 



1 

27 



2-7T [ m 



1 



± 



(18) 



the sign in front of the phase depending on the direc- 
tion of the supercurrent flow. Here r = vp/L is the 
time of flight between two successive Andreev reflections 
and L is the corresponding length of the trajectory. The 
supercurrent-carrying density of states is then found from 



js(E;< 



oc 



(19) 



resulting into a peaklike spectrum that contains states 
carrying both positive (E = E+) and negative (E = E^) 
supercurrent. In the presence of disorder, the distribu- 
tion of the times of flight r depends on the impurity 
potential and the spectral supercurrent is conveniently 
characterized by its impurity-averaged smooth density of 



states. However, the resulting spectrum still contains 
many properties similar to the clean limit, such as the 
varying sign of the supercurrent carried at different ener- 
gies. This analogy holds, even though on a formal level, 
the calculation within our quasiclassical technique does 
not directly invoke these concepts. 



III. SHORT- AND LONG-JUNCTION LIMITS 

The spectrum of current-carrying states in the weak 
link depends very much on the ratio of the length L of 
the wea k link and the superconducting coherence length 
£o = y/D/2A, or in other words, on the ratio between 
the superconducting order parameter A and the Thou- 
less energy Et = D/L 2 of the weak link. In the case of a 
long weak link, 1 L>^ ( or 7 equivalently, Et <C A), the 
spectrum is wide and many energy states contribute to 
the supercurrent with only a small phase-dependent gap 
of the order of a few Et at low energies. In the opposite 
limit, only the states with energy E G [A| cos(^/2)|, A] 
carry supercurrent and between these limits, A serves 
as a cutoff for the spectral supercurrent: there may ex- 
ist some current-carrying states with E > A, but their 
contribution vanishes quickly with E — A. 

The energy ranges can be understood as follows: In the 
center of the junction, which is the bottleneck for the su- 
percurrent, both superconductors provide sufficient cor- 
relations that a gap in the energy spectrum of a size E g 
is induced at <j> — 0, where the size of E g interpolates 
between Et (long junction) and A (short junction). 31 ' 32 
If now a finite phase difference is applied, the correla- 
tions from either side start to interfere more and more 
destructively leading to a closing of the gap at <f> = 7r. 33 
Hence the lower energy bound, below which no bound 
states exist, is set by this phase-dependent gap. Above 
A, the states depend less and less on the superconduct- 
ing properties, hence their phase dependence is rapidly 
lost and they also do not contribute to the supercurrent. 



A. Short-junction limit L <C £o 

In the limit when the superconducting order parameter 
A is much smaller than the Thouless energy, the super- 
current is carried by states with energies much below Et 
and we may thus neglect the first term on the right side of 
Eq. (5). In this case, we get an analytical solution to the 
differential equations without further approximations, 



0(x) 



( Va 2 + 1 \ 
arcosh cosh[j E a(x — x )] , (20) 



x( x ) — Xo — arctan{atanh[jsa(x — #o)]}> 



(21) 



where a and Xq are constants which along with the spec- 
tral supercurrent jp are determined from the boundary 
conditions. In the two-probe case we can choose the 
origin in the center of the weak link, and assume the 



■5 



functions 9(x) and x( x ) take the bulk values at the NS 
boundary (x = ±L/2). Thus we get x a = and 



Je 



y/E 2 - A 2 cos 2 (0/2) 
Acos(0/2) : 

2Acos(</>/2) 
^JE 2 - A 2 cos 2 (</>/2) 



(22) 



arcosh 



/£ 2 - A 2 cos 2 (</>/2) 



£ 2 - A 2 



(23) 



In the real-time calculation of the supercurrent, we are 
mostly interested in the imaginary part of the spectral 
supercurrent ]e- This is 



'0, E>A 

7rAcos(0/2) 
y/E 2 -A 2 cos 2 (0/2 

,0, |£7| < A|cos(0/2)| 



Ee [A|cos(<£/2)|,A] 



(24) 

and Im[j E (-E)} = -lm[j E {E)]. At T = 0, we get the ob- 
servable supercurrent by simply integrating Im(js) over 
the energy to obtain 



7rA cos(</>/2) , , . , , 
7 s = _ n ; artanh[sin((/)/2)]. 



eR 



N 



(25) 



For a finite temperature, we have to multiply this by 
the distribution function tanh(_E/2fc^r) and integrate 
over the energy, which is conveniently done using the 
Matsubara technique [i.e., substituting E = iiv n , where 
uj n = nT(2n + 1) are the poles of tanh(.E/2fc B T) and 
summing RcQ'b) over n = 0, 1, see Refs. 10 and 20 for 
details], yielding 



2ttAT 
Is = — - — cos 



eR 



N 



x arctan 



1 



^ VA 2 cos 2 (0/2)+c 2 
Asin(</>/2) \ 



(26) 



^A 2 cos 2 (4>/2)+uj 2 I ' 



As expected, Eqs. (25) and (26) are the same as obtained 
by Kulik and Omel'yanchuk 34 and in numerical studies 10 
for the same limit. 

In a setup where the distribution function can be con- 
trolled by an additional probe coupled to the system via a 
narrow normal wire (such that the current-carrying states 
are not essentially deformed), the resulting supercurrent 
as a function of control voltage V at T = reads 



ttAcos(</>/2) 
2eR N 



In 



A(l + sin(</>/2)) 



V+^/V 2 - A 2 cos 2 (0/2) 

(27) 

for V <G [A|cos(0/2)|, A]. Above A, Is vanishes, and for 
V < A|cos(</>)|, the supercurrent has the form of Eq. (25), 
independent of V. 

The spectral supercurrent of Eq. (24) can also be ob- 
tained from the diffusive limit of the corresponding quan- 
tity derived in Ref. 35. There, the supercurrent is written 




FIG. 2: Spectral supercurrent for a few values of A/Et- (a) 
Junctions longer than the superconducting coherence length 
£o- A finite A/Et shows up as a peak at E = A. (b) 
Short junctions L ~ £o: the peak at E = A persists, but 
another develops around E — Acos(0/2). In (b), the spec- 
tral supercurrent is normalized by Et/A to allow for the an- 
alytical solution at Et/A — > oo. The inset shows how the 
zero-temperature, zero-voltage critical current behaves as a 
function of A/Et, in accordance with Ref. 10. 



as a sum of the contributions from different bound states, 

N 

■ sii i ( < o 1 



/ s = ^ sinW ^^tanh^ 
p=i p 



2k B T J ' 



(28) 



where the bound-state energies E p depend on the trans- 
mission eigenvalues r p by E p = A[l — t p sin 2 (0/2)] 1 / 2 . 
Writing Eq. (28) in the form of an energy integral, 

I s = e A sH(j)) JdEf2%^(-^f)s(E-E p ), 

(29) 

and averaging the transmission eigenvalues 
over their diffusive-limit distribution, 36 p(r) = 
(7r/2e 2 i?Ar)(rv / l — t) -1 , yields a spectral supercur- 
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FIG. 3: Energy dependence of the Fourier sin-transformed 
spectral supercurrent: four lowest harmonics n = 1,2,3,4 
corresponding to the phase dependencies sin(n0). The energy 
scales and magnitudes of the different harmonic constituents 
decay as 1/n. 



rent given by Eq. (24) multiplying the distribution 
function ta,nh(E/2k B T). 
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FIG. 4: Voltage dependence of the amplitudes of the first four 
harmonics of the observable supercurrent. Inset: Correspond- 
ing temperature dependence (even harmonics yield negative 
amplitudes, but here we plot the absolute values of the super- 
current). 



IV. PHASE DEPENDENCE 



B. Long and intermediate-length junctions 

If the length L of the weak link is much longer than 
£o, the supercurrent is carried by a wide spectrum of 
energies 1,2,37 . At low E, however, the current-carrying 
density of states has a phase-dependent minigap rem- 
iniscent of the gap in the usual density of states of a 
SNS sample 29 . Above the gap, ImQ'^) rises sharply, 
then starts to oscillate with an exponentially decay- 
ing envelope. This oscillatory behavior is responsible 
for the occurrence of the tt state in nonequilibrium- 
controlled Josephson junctions 5,6 and in ferromagnetic 
weak links. 38,39 

The spectral supercurrent as a function of energy is 
plotted for a few values of A/Et in Fig. 2, the upper 
figure showing the limit Et <C A and the lower the limit 
E T ^ A. 

For a finite ratio A/Et, the divergence of the density 
of states at the superconducting gap edge is reflected as 
a peak in the spectral supercurrent at E = A. The di- 
rection of the peak, positive or negative, is determined 
by geometric considerations and hence depends on the 
precise value of A/Et- For Et ~ A, the spectral super- 
current \m{jE) tends towards the short-junction result, 
Eq. (24), replacing the Thouless gap by the gap of width 
Acos(<?!>/2). Moreover, for A/E T -> 0, the width of the 
peak at E = A tends to zero. 

In the limit T » Et for a long junction (Et <C A), 
the temperature dependence of the obtained observable 
supercurrent tends to the limits considered by Likharev 40 
for A <C T and by Zaikin and Zharkov 41 for a general 
A/T. 



Originally, the Josephson effect was discovered for in- 
sulating weak links in the tunneling regime, i.e., in the 
limit of a very low tunneling probability. There, the su- 
percurrent is due to an uncorrelated transfer of Cooper 
pairs through the weak link. 42 As a result, one obtains 
the familiar dc Josephson relation Is = I c sm((j)). How- 
ever, it has been shown (see, e.g., Refs. 13 and 34) that 
other kinds of weak links, through which the transmis- 
sion probability is much above zero, may have a different 
current-phase relation. Thus we may write in general 

oo 

7c7(0) = ^/5sin(n0), (30) 

n=l 

where the amplitudes are the coefficients of the 
Fourier sine series of /c(<^)- Fot example, in a ballistic 
weak link where the transmission probability for Cooper 
pairs is 1, Iq oc — (— l)™/n, yielding a sawtooth form for 
Jc(</>)- 13,16 The odd parity with respect to <f) (appearance 
of only sine terms) of this representation reflects the fact 
that the supercurrent is driven by the spatial asymmetry 
introduced by the application of the phase: changing the 
phase to a negative value corresponds to mirroring the 
structure about the center and hence to a reversal of the 
current. 

The occurrence of higher harmonics in Eq. (30) may 
be interpreted as a correlated transfer of n Cooper pairs 
through the weak link as a result of the pairing corre- 
lations extending through quasiparticle paths containing 
multiple Andreev reflections. For example, the flux quan- 
tum for the nth harmonic is h/2ne, i.e., corresponding to 
a charge 2ne. 

In a diffusive weak link considered in this paper, the 
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transmission probabilities for the Cooper pairs are widely 
distributed between zero and l. 36 As a result, one may 
get contributions from the higher harmonics to the phase 
dependence. This is shown in Fig. 3 where the ampli- 
tudes Im(j||) of the first four harmonics of the Fourier 
sine transformed spectral supercurrent through a long 
weak link (L 3> £o) are plotted as a function of en- 
ergy. Both the amplitude of Im(f^) and the effective 
energy scales decay with a power of n suggesting that 
the observation of the higher harmonics is easiest at low 
temperatures. The corresponding temperature and volt- 
age dependencies of the critical currents plotted in 
Fig. 4 behave analogously. A numerical fit to the ob- 
tained critical currents at eV = ksT = yields roughly 
ig oc -{-l) n /n 2 and to the voltage V* where Ig(V) first 
changes sign suggests that the effective energy scales be- 
have as E* = Et(ci + Cijn), with some constants c\p,. 

This behavior can be understood by identifying the 
higher harmonics with the correlated transfer of a clus- 
ter of n Cooper pairs. Now instead of the phase <f>, the 
cluster has the phase n<p and since the cycle contains 
2n Andreev reflections, the effective trajectory length is 
increased from L to nL. In Eq. (18), making these re- 
placements yields the observed result, E* oc (ci + cijn). 
In the diffusive limit, the effective trajectory length in- 
creases in the second power of the length of the weak link, 
l e H oc L 2 , but since the phase is reset after each traver- 
sal through the weak link, we simply get Z e ff,n oc nL 2 . 
Hence, similarly to the alternating sign, the scaling of 
the effective energies with index n follows the behavior 
of the ballistic- limit spectral supercurrent. 

Since the crossover voltages V*, where /g(V^) = 0, 
depend on n, the actual critical current never vanishes 
at the crossover: it is rather that the current-phase rela- 
tion changes its form near the crossover voltages. Such a 
change was observed in Ref. 11, where the current-phase 
relation of a controllable Josephson junction was mea- 
sured in a superconducting quantum interference device 
(SQUID) geometry. 

In the short-junction regime L <C £o the contributions 
of the different harmonics can be derived analytically. A 
general form for Im(jjj) would be complicated, but as an 
example, the first two amplitudes are 



Ini(.?i) = 
Mi!) = 



E 
A 

£ 2 (2A 2 - 3£ 2 ) 
A 3 ' 



(31) 
(32) 



Clearly all the harmonics share the same energy scale, 
A, but as for the case of a long weak link, the amplitude 
decays also here roughly as 1/n 2 . Namely, integrating 
Im(jg) over the energy, we get for a short weak link 



(-l) n eA 



R N {2n+ l)(2n- 1) 



(33) 




Replacing A by Et and scaling by a numerical factor 
close to 33, this form fits the amplitudes of the long- 



FIG. 5: Spectral supercurrent for different cross sections A c 
of the control probe. From top to bottom, the cross sec- 
tion is A c — 0, 0.2, 1, 2, 4 times the cross section of the weak 
link. Inset: observable supercurrent as a function of A c for 
T = V = (upper set of circles) and for ksT = 3Et 
(lower set of circles). The zero-temperature result is fitted 
to I,s(Ac = 0)A w /(A w + Ac/3) (solid line) and the finite- 
temperature result to Is(A c — 0)A w /(A w + A c /2) (dashed 
line). 



junction harmonics as well but a rigorous proof does not 
exist. 



V. EXTRA TERMINALS 

In order to relate our results to physical observables, 
we have to evaluate statistical expectation values. In two- 
probe SNS weak links in equilibrium, most of the experi- 
mental observations have been accurately described with 
the equilibrium Matsubara technique. 10 However, one of 
the recent advances in the research of the Josephson ef- 
fect has been done in nonequilibrium situations where the 
distribution function in the weak link has been controlled 
by coupling one or more normal-metal reservoirs to the 
weak link by phase-coherent wires. 5,6 While making it 
possible to control the occupation of the current-carrying 
states, these extra wires also affect the form of Im(js)- 2 
In the discussion of these effects, we concentrate on the 
regime of a long junction, L 3> Co- There, most notably, 
the control probes allow for the existence of states with 
low energies, and therefore the Thouless gap is lifted. 
Moreover, the existence of the normal reservoirs brings 
some extra depairing by imposing a vanishing boundary 
condition for the pairing amplitude / = sinh(0) inside 
the reservoirs. As a result, the amplitude of the spec- 
tral supercurrent decreases due to the extra probes. In 
what follows, we consider the effects of the integrated 
cross-sectional area A c of the control wires attached to 
the center of the weak link with cross section A w (note: a 
similar effect would be present if the control wires and the 
weak link were made of different materials with normal- 
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4 5 10 15 20 



eV/E T 

FIG. 6: Voltage dependence of supercurrent for <j> — it /2 in 
the presence of the control probe with different cross sections 
A c . Inset: The corresponding temperature dependence of the 
supercurrent. The cross sections have been chosen as in Fig. 
5. 

state conductivities <jn,c and a n, w — however, here we 
simply talk about A c and A w ) and of the length L c of 
the phase coherent control wires, compared to the length 
L of the weak link. For simplicity, we assume that the 
widths of the control wires are much smaller than L, al- 
lowing one to treat these wires as quasi-lD structures, 
connected by the rules of Nazarov's circuit theory 25,36 . 
In the language of this circuit theory, the extra normal 
wires divert some of the spectral current to the normal 
reservoirs, thus decreasing the pairing correlations be- 
tween the superconductors. 

The spectral supercurrent as a function of energy for 
different cross sections A c of the control wire is plotted 
in Fig. 5. Here we have taken the length L c — hL 3> L. 
Already a small A c <C A w yields a finite Im(j£) at low 
energies, but does not much reduce the total magnitude 
of the supercurrent. The resulting temperature and volt- 
age dependencies of the total supercurrent I s is plotted in 
Fig. 6. Except at the very lowest voltages/temperatures 
of the order of the Thouless gap, the extra probes do not 
change the voltage/temperature dependence of I s from 
the two-probe case, but only the overall magnitude is 
decreased. From the resulting Is{A c ) we obtain 

I S (A W ,A C ) = — ^ W A ^ 2 I S (A W ,A C = 0), (34) 

which holds very well for max(eV, ksT) ^> Et- 

If the length L c of the phase-coherent control wire is 
larger than L, the effect of the control wires is indepen- 
dent of the precise value of L c . For L c < L, the spec- 
tral supercurrent is altered for E < TiD/L"^ such that 
the overall magnitude is decreased, the observable super- 
current tending towards zero as L c ^ 0. The spectral 
supercurrent and the resulting temperature/voltage de- 
pendencies for A c = A w and for three different L c /L are 




FIG. 7: Spectral supercurrent for different lengths L c of the 
control probe. From top to bottom, the length is L c = 
5,0.25,0.15 times the length of the weak link. The cross 
section of the control probe is chosen equal to that of the 
weak link. Inset: Supercurrent at T = (solid line) and 
T — 3Er/kB (dashed) as a function of the length L c of the 
control probe with respect to the length of the weak link. 
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4 5 10 15 20 

eV/E T 

FIG. 8: Voltage dependence of supercurrent for cj> — ir/2 in 
the presence of a control probe with length L c and cross sec- 
tion equal to that of the weak link. Inset: The corresponding 
temperature dependence of the supercurrent. The lengths L c 
have been chosen as in Fig. 7. 



plotted in Figs. 7 and 8, respectively. 



VI. NONIDEAL INTERFACES 

Normal-metal — superconductor interfaces with re- 
duced transmissivity can be taken into account by spec- 
ifying the transmission eigenvalues through the interface 
and taking them into account as in Eq. (7). Here we 
consider a typical case described by the distribution of 
eigenvalues for a "dirty" interface. 
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0,5 5 10 15 20 25 30 35 40 
E/E T 

FIG. 9: Spectral supercurrent in a weak link with dirty in- 
terfaces to the superconductors, characterized by the ratio 
Tb = Ri/Rn of the resistances. The interface resistance 
Ri is assumed the same for both interfaces. From top to 
bottom: r b = 0, 0.2, 0.4, 0.6, 0.8, 1.0. Enhanced scatter- 
ing at the interface reduces mostly the amplitude of the su- 
percurrent, but also slightly the effective energy scale. In- 
set: zero-temperature, zero-voltage supercurrent as a func- 
tion of r b . Circles: calculated supercurrent; solid line: fit to 
Is{n) = Is(r b = 0)/(l + 1.6r b ). 




5 10 15 20 



eV/E T 

FIG. 10: Voltage dependence of the supercurrent with dirty 
NS interfaces at <f> = 7r/2. Inset: Corresponding temperature 
dependence. The values of Ri are the same as in Fig. 9. 

Figure 9 shows the spectral supercurrent for a long 
junction connected to superconductors through a dirty 
interface with resistance Ri (yielding a total resistance 
Rn + 2Rj between the superconductors). Due to the 
additional resistance, the amplitude of the supercurrent 
decays with rb, but also the energy scales decrease since 
the interface barrier can to certain extent be thought as 
adding a barrier-equivalent length 43 to the path length of 
the quasiparticles. Observing the temperature and volt- 
age dependencies of the resulting supercurrent, plotted 
in Fig. 10, shows that the amplitude of the supercurrent 



behaves at T = 0, V = as 

Ic(Ri,Rn)~ Rn + 16Ri - (35) 

i.e., the resistances should not simply be added up but 
the dirty interface decreases the supercurrent less effi- 
ciently than the normal-metal resistance. Furthermore, 
the effective energy scale E* found, e.g., from the voltage 
dependence indicates that it follows the approximate law 
E* = E T (1 + OJRj/Rn). This kind of a behavior of the 
critical current and the spectral supercurrent is similar 
to those found in Ref. 24 [especially, see Eqs. (55) and 
(56)] and Ref. 2 in the case of a tunneling interface. 

VII. CONCLUSIONS 

In this paper, we have systematically investigated the 
spectrum of current-carrying states in a phase coherent 
normal-metal weak link. Taking into account the effect 
of extra terminals, the characteristic energy scales — the 
Thouless energy Ex and the BCS superconducting gap 
A — and a finite NS interface resistance makes it easier 
to find a quantitative agreement with the obtained ex- 
perimental results on the noncquilibrium-controlled su- 
percurrent. We have also been able to derive analytical 
results in a number of limits. Moreover, we have dis- 
cussed the underlying microscopic phenomena leading to 
the 7T state and have explained its properties, such as its 
dependence on energy, and higher harmonics in the phase 
dependence, by invoking Andreev bound states smeared 
over a broad distribution of times of flight, and by mul- 
tiple Andreev cycles tranferring more than one Cooper 
pair in a single coherent process. 

To obtain an optimal voltage control of the supercur- 
rent, the interface resistances should be much smaller 
than the weak-link wire resistance, the control wire 
should be slightly longer than the weak link (but not 
much longer to reduce the inelastic effects on the nonequi- 
librium distribution function), and as thin as possible 
compared to the weak link. For the observation of the 7r 
state, the ideal limit is the long-junction limit L > ^o, 
where the spectrum of the current-carrying states is not 
cut off by the superconducting gap. 

In typical experiments, also inelastic scattering ne- 
glected here may become important. However, since 
equilibrium phenomena induced by the superconducting 
proximity effect have been quantitatively described by 
the quasiclassical theory without incorporating such ef- 
fects (see, e.g., Rcfs. 10 and 44), we expect these inelastic 
effects to be mostly important in the kinetic equations 
describing the noncquilibrium distribution functions. In 
recent experiments, 6 ' 8 including these inelastic terms into 
the kinetic equations has lead to good agreement between 
the theory and the experiments. Therefore our results 
provide an accurate and independent way of also char- 
acterizing such inelastic effects by observing how they 
affect the supercurrent. 
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